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Critical mass of bacterial populations in a generalized Keller-Segel model. 
Analogy with the Chandrasekhar limiting mass of white dwarf stars 
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We point out a remarkable analogy between the limiting mass of relativistic white dwarf stars 
(Chandrasekhar's limit) and the critical mass of bacterial populations in a generalized Keller-Segel 
model of chemotaxis [Chavanis & Sire, PRE, 69, 016116 (2004)]. This model is based on generalized 
stochastic processes leading to the Tsallis statistics. The equilibrium states correspond to polytropic 
configurations similar to gaseous polytropes in astrophysics. For the critical index — d/(d — 2) 
(where d > 2 is the dimension of space), the theory of polytropes leads to a unique value of the 
mass M c that we interpret as a limiting mass. In d — 3, we find M c = 202.8956... and in d = 2, we 
recover the well-known result M c = 8n (in suitable units). For M < M c , the system evaporates (in 
an infinite domain) or tends to an equilibrium state (for box-confined configurations). For M > M c , 
the system collapses and forms a Dirac peak containing a mass M c surrounded by a halo. This 
paper exposes the model and shows, by simple considerations, the origin of the critical mass. A 
detailed description of the critical dynamics of the generalized Keller-Segel model will be given in a 
forthcoming paper. 

Keywords: Chemotaxis; Generalized thermodynamics; Nonlinear meanfield Fokker-Planck equations; Self- 
gravitating Brownian particles 



I. INTRODUCTION 

Recently, there has been a growing interest for the dy- 
namics and thermodynamics of systems with long-range 
interactions Tj. Such systems are numerous in Nature 
and share fascinating analogies. For example, the statis- 
tical mechanics of stellar systems (like globular clusters 
and elliptical galaxies) in astrophysics and the statistical 
mechanics of large-scale vortices (like Jupiter's great red 
spot) in two-dimensional turbulence present deep simi- 
larities despite the very different physical nature of 
these systems. Some connections have also been devel- 
oped between the dynamics of stellar systems and the 
dynamics of the Hamiltonian Mean Field (HMF) model 
[3, lij and of the free electron laser [f| . More recently, the 
authors have studied a model of self-gravitating Brow- 
nian particles [f| 0, HI and discovered striking analogies 
with the Bose-Einstein condensation in the canonical en- 
semble [91 and with the phenomenon of chemotaxis in 
biology [1 01 ] . In this paper, we point out a novel analogy 
between the Chandrasekhar limiting mass of white dwarf 
stars ll| and the critical mass of bacterial populations 



experiencing chemotactic aggregation. 

The name chemotaxis refers to the motion of organisms 
induced by chemical signals [12]. In some cases, the bio- 
logical organisms (bacteria, amoebae, endothelial cells, 
ants...) secrete a substance (pheromone, smell, food, 
...) that has a long-range attractive effect on the organ- 
isms themselves. Therefore, in addition to their diffu- 
sive motion, they move preferentially along the gradient 
of concentration of the chemical they secrete (chemotac- 
tic flux). When attraction prevails over diffusion, the 
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chemotaxis can trigger a self-accelerating process until a 
point at which aggregation takes place. This is the case 
for the slime mold Dictyostelium discoideum and for the 
bacteria Escherichia coli. This is referred to as chemo- 
tactic collapse. A model of slime mold aggregation has 
been introduced by Keller & Segel [IH in the form of 
two coupled differential equations. A simplified version 
of this model has been extensively studied in the case 
where the degradation of the secreted chemical can be 
neglected and an immediate production is assumed [l4[. 
In that case, the Keller-Segel (KS) model become isomor- 
phic to the Smoluchowski-Poisson (SP) system describing 
self-gravitating Brownian particles [tJ . The steady states 
correspond to isothermal distributions similar to isother- 
mal stars in astrophysics [l5[. On the other hand, the 
KS model and the SP system can be viewed as standard 
mean field Fokker-Planck equations associated with the 
Boltzmann free energy [l6|, [T3] ■ In this sense, they are 
based on ordinary thermodynamics. Since these systems 
are strongly dissipative, the correct statistical ensemble 
is the canonical ensemble Q- 

Recently, modified forms of Keller-Segel models have 
been introduced in order to describe more general sit- 
uationsjll E3, H3, HE HI, HE E3|- For example, in 
Ref. [20J, we have introduced and studied a general- 
ized Keller-Segel (GKS) model of chemotaxis taking into 
account anomalous diffusion. In this model, the coeffi- 
cient of diffusion is assumed to depend on the density 
like D(p) = Dp 1 ^ 1 with 7 = 1 + 1/n (the ordinary KS 
model is recovered for 7 — 1 or n — > +00). Anomalous 
diffusion is known to appear in many problems of biology 
[l2| and it is likely that it can play a role in the process of 
chemotaxis. The generalized Keller-Segel (GKS) model 
is isomorphic to the generalized Smoluchowski-Poisson 
(GSP) system describing self-gravitating Langevin par- 
ticles [2fJ. The steady states correspond to polytropic 
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distributions similar to polytropic stars in astrophysics 
On the other hand, the GKS model and the GSP 
system can be viewed as nonlinear mean field Fokker- 
Planck equations (2f| [26[ associated with the Tsallis free 
energy [23 ■ In this sense, they are related to a notion of 
(effective) generalized thermodynamics. 

For the GKS model, there exists a particular polytropic 
index n.3 = d/{d— 2) (where d > 2 is the dimension of 
space) at which the dynamics is critical. For this index, 
the theory of polytropes leads to a unique value of the 
mass M c , independent on the size of the object, that 
we interpret as a limiting mass. This is the counter- 
part of the Chandrasekhar limiting mass [TlJ for white 
dwarf stars which are equivalent, in the ultra-relativistic 
limit, to polytropes of index n — 3 (in d = 3). This 
unexpected analogy between two systems that have ap- 
parently nothing in common was pointed out i n 12811 and 
is here systematically developed (see also Ref. |29||). 

In this paper, we expose the GKS model (using the 
notations of biology) and show its relation to Tsallis gen- 
eralized thermodynamics. Then, we show by simple con- 
siderations, the origin of the critical mass M c appearing 
at the polytropic index 713 = d/(d — 2). For d = 3, we 
find that M c = 202.89561... in a proper system of units. 
For d — 2 where 713 — > +00, we recover the critical mass 
M c = 8tt corresponding to the ordinary KS model in two 
dimensions whose critical dynamics has been extensively 
studied in 0, 0, El (see (3fl, M, H [H, H HE H for 
many rigorous results obtained by applied mathemati- 
cians). For d = 3, a detailed description of the critical 
dynamics of the GKS model for M < M c and M > M c 
in bounded and unbounded domains will be given in a 
forthcoming paper [29(. 



II. THE ORDINARY KELLER-SEGEL MODEL: 
BOLTZMANN THERMODYNAMICS 

The primitive Keller-Segel model KUp describing the 
chemotaxis of bacterial populations [56j consists in two 
coupled differential equations 



dp 



= V • (D 2 Vp) - V • (£>iVc) 



dc 
di 



D c Ac+ h{c)p~ k(c)c, 



(1) 



(2) 



which govern the evolution of the density of bacteria 
p(r,t) and the evolution of the secreted chemical c(r,t). 
The bacteria diffuse with a diffusion coefficient D2 and 
they also move in a direction of a positive gradient of 
the chemical (chemotactic drift). The coefficient D\ is a 
measure of the strength of the influence of the chemical 
gradient on the flow of bacteria. On the other hand, the 
chemical is produced by the bacteria with a rate h(c) and 
is degraded with a rate k(c). It also diffuses with a dif- 
fusion coefficient D c . In the general Keller-Segel model, 
D\ = Z?i(p, c) and D2 = -^(p, c) can both depend on 



the concentration of the bacteria and of the chemical. 
This can take into account microscopic constraints, like 
close-packing effects or anomalous diffusion [iij . 

A simplified version of the Keller-Segel model is pro- 
vided by the system of equations 



dt 



(DVp - xpVc) 



dc 
di 



D c Ac + hp — kc, 



(3) 



(4) 



where the parameters are positive constants. Another 
simplification is obtained in a limit of large diffusivity of 
the chemical D c — > +00 and for sufficiently large concen- 
trations of the bacteria (see [l4| and Appendix C of 24 1 
for details). In that case, Eq. ^ is replaced by a Poisson 
equation 



-Xp, 



(5) 



where A = h/D c . Equation ^ can be viewed as a mean- 
field Fokker-Planck equation associated with a Langevin 
dynamics of the form 



rfr 



= X Vc+V2DTL(t), 



(6) 



where R(i) is a white noise, D is a diffusion coefficient 
and x plays the role of a mobility. The Langevin equation 
describes a point organism performing a random walk 
biased in the direction of a drift velocity proportional to 
the local chemical gradient. 

The KS model ©-([I]) monotonically decreases the 
Lyapunov functional 



F 



pc dr 



D 



p In p dr, 



(7) 



which is similar to a free energy F — E — T e fjS where 
E = —(1/2) J pcdr is the energy of interaction and S — 
— J phipdr is the Boltzmann entropic functional. We 
have defined an effective temperature 



T - D 

1 eff - — 1 



X 



(8) 



which is given by a form of Einstein's formula. Using the 
KS model, we find that 



1 

XP 



{DVp-xP^cfdv < 0, 



(9) 



which is similar to the proper version of the H-theorem in 
the canonical ensemble [13, HE] • The stationary solutions 
of Eq. ([3]), corresponding to F = 0, are given by 



p = Ae° c , 



(10) 



where A is determined by the mass M. This is similar 
to the Boltzmann distribution for a system in a potential 
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— c(r) at temperature T e /y = Dj\. The steady states of 
Eq. are critical points of the free energy F at fixed 
mass M. They cancel the first order variations 6F — 
aSM — where a is a Lagrange multiplier. Moreover, 
it can be shown [25| that they are linearly dynamically 
stable for the KS model if and only if they are (local) 
minima of F at fixed mass. The equilibrium state is 
obtained by coupling Eq. (|10p to Eq. JS]) leading to the 
Boltzmann-Poisson equation. 

These analogies with thermodynamics take even more 
sense if we remark that the Keller-Segel model is isomor- 
phic to the Smoluchowski-Poisson system 



dt 



k R T. 



-Vp + pV$ 



A$ = S d G P , 



(11) 



(12) 



describing the dynamics of self-gravitating Brownian par- 
ticles in an overdamped limit and in a mean field approx- 
imation 0, ■ We have the correspondence 

D = k B T/tm, x = Ve, c = -$, \ = S d G. (13) 

We note, in particular, that the role of the gravitational 
potential $ in the SP system is played, in the KS model, 
by the concentration — c of the chemical. Therefore, the 
pheromonal substance has a long-range attractive effect 
similar to the gravitational attraction in astrophysics. In 
biology, the chemotactic aggregation is mediated by a 
physical field (the chemical produced by the organisms) 
while the nature of the gravitational force in astrophysics 
is more abstract. The equilibrium states of the SP system 
correspond to a condition of hydrostatic balance 



VP + pV$ = 0, 



A$ = S d G P , 



with an isothermal equation of state 
P = pkBT/m. 



(14) 



(15) 



Therefore, the equilibrium states of the KS model and 
SP system have the same structure as isothermal stars 
in astrophysics [l5| . Of course, the dynamics of the KS 
model and SP system is different from the dynamics of 
stars which is rather described by hydrodynamic equa- 
tions like the Euler-Poisson system [33] . 



III. THE GENERALIZED KELLER-SEGEL 
MODEL: TSALLIS THERMODYNAMICS 

In a previous paper [2(| (see also [Hj]), we have in- 
troduced and studied a generalized Keller-Segel model of 
the form 



dp 
dt 



V • (DVp-> - xpVc) , 



Ac = -Xp. 



(16) 



(17) 



Comparing with the primitive Keller-Segel model JU-JU, 
it corresponds to D% = D^p 1-1 and D\ = \P- For sim- 
plicity, we have considered that the chemical is deter- 
mined by a Poisson equation (|17p . More generally, we 
can consider that it satisfies the field equation J4|) . How- 
ever, the analytical results that we obtain in Sees. IIVI 
and [V] are only valid for the Poisson equation Q17[) . For 
7=1, we recover the ordinary Keller-Segel model ([5])- 
J5]). Equation ([16]) can be viewed as a nonlinear mean 
field Fokker-Planck equation of the form considered in 
[25l [26l [H, [39[ . It can be obtained from the stochastic 
process 



dr 
~dt 



(18) 



where R(£) is a white noise [57j . This equation describes 
a situation where the mobility x is constant but the dif- 
fusion coefficient D(p) = Dp 1-1 can depend on the den- 
sity. This can account for anomalous diffusion and non- 
ergodic behaviors. For 7=1, we recover the ordinary 
Langevin equation J6j with constant diffusion coefficient 
D and constant mobility x- The stochastic process (fT8|) 
has been introduced by Borland 14211 in relation with Tsal- 
lis generalized thermodynamics [271 ]. For 7 = 1, we have 
a pure random walk. In that case the sizes of the ran- 
dom kicks are uniform and do not depend on where the 
particle happens to be. For 7^1. the size of the ran- 
dom kicks changes, depending on the distribution of the 
particles around the "test" particle. A particle which is 
in a region that is weakly populated [small p(r, t)] will 
tend to have smaller kicks if 7 > 1 and larger kicks if 
7 < 1. Since the microscopies depends on the actual 
density, this creates a bias in the ergodic behavior of the 
system. 

The GKS model decreases the Lyapunov functional 



If, D 1 
F = — / pedv H 

iy X7-1 



(p 7 - p)dr. 



(19) 



It can be interpreted as a generalized free energy of the 
form F = E-T e ffS where E = -(1/2) / pedv is the en- 
ergy and S = —1/(7—1) /(p 7 — p)dr is the Tsallis entropy 
(the polytropic index 7 plays the role of the Tsallis q pa- 
rameter) . The effective temperature T e f f is still given by 
the Einstein- like formula J8|). Using the GKS model, we 
find that 



F = 



1 

XP 



(PVp 7 - xP^cfdv < 0, 



(20) 



which is similar to the proper version of the ii-theorem in 
the canonical ensemble in a generalized thermodynamical 
framework [l7], HE] • The stationary solution of Eq. JT| 
corresponding to F = 0, is given by 



P : 



fi 



X 7- 1 



D 



7 



1/(7-1) 



(21) 



where p is determined by the mass M. This corresponds 
to the Tsallis distribution with inverse temperature (3 = 



4 



1/T e ff = x/D and "g-parameter" 7. The steady states 
of Eq. (|16[) are critical points of the free energy F at 
fixed mass M. They cancel the first order variations SF — 
aSM — where a is a Lagrange multiplier. Moreover, it 
can be shown that they are linearly dynamically stable 
for the GKS model if and only if they are (local) minima 
of F at fixed mass (2f| . The equilibrium state is obtained 
by coupling Eq. (JUJ) to Eq. (JT7J) leading to the "Tsallis- 
Poisson" equation. 

These analogies with generalized thermodynamics take 
even more sense if we remark that the GKS model is iso- 
morphic to the generalized Smoluchowski-Poisson system 



dt 



~ [K\Jp~< + pV$) 



(22) 



A$ = S d G P , (23) 

provided that we set 

D = K/(, x = Vf, c = -$, X = S d G. (24) 

In particular, the equilibrium states correspond to a con- 
dition of hydrostatic balance 

VP + pV$ = 0, A$ = S d Gp, (25) 

with a polytropic equation of state 

P = Kp~ f . (26) 

Therefore, the equilibrium states of the GKS model and 
GSP system have the same structure as polytropic stars 
in astrophysics fl5T| . As in astrophysics, we define the 
polytropic index n by 



1 

7 = 1 + -. 

n 



(27) 



where po is the central density, and using the Poisson 
equation (jTTJ) with the steady distribution (|2"Tj) , we find 
after simple algebra that the function 0(£) is solution of 
the Lane-Emdcn equation 



1 



(29) 



with 6 — 1 and 9' = at £ = 0. In the following, we 
consider d > 2 and n > 0. It is shown in [2(| that 
polytropic configurations are self-confined iff 



n < n 5 



d + 2 
d-2' 



(30) 



In that case, the function vanishes at a finite distance 
£ = £1. Consequently, the density p(r) vanishes at a 
distance i?* defining the radius i?* of the polytrope. The 
relation between the radius and the central density is 



\ 1-1/r. 
D(l + n) 



1/2 



R*. 



(31) 



The mass M — J^* pSd,r d 1 dr of the configuration is 
given by 



M = S dP o 



' D(l + n)~ 


d/2 


\ 1-1/n 


/ 



(32) 



Using the Lanc-Emden equation (|29)) . we get 

ci/2 



M 



-SdPo 



D(l + n) 

T 1-1/n 



(33) 



Expressing the central density in terms of the radius, 
using Eq. (|3ip . and introducing the index 



IV. THE CRITICAL MASS OF BACTERIAL 
POPULATIONS FOR n = n 3 

The steady states of the GKS model are determined 
by substituting the polytropic distribution (|2ip in the 
Poisson equation pT|) . The resulting configurations have 
the same structure as polytropic stars in astrophysics. 
Therefore, we can readily apply the theory of polytropes 
developed long ago by Emden [43| to the biological prob- 
lem. For sake of generality, we extend these results to a 
space of d dimensions [2(J. The dimension d — 2 (which 
is important in biology) will be considered specifically in 
Sec. El 

We consider spherically symmetric configurations. 
Defining 



-l/n 



1 1/2 



D(l + n) 



r, 



(28) 



n 3 



d-2' 

we obtain the mass-radius relation 



(34) 



M (n-l)/n R [(d-2)(n 3 - n)] /n = gft+g) ^(n-D/n (35) 

Ax ^(l-n)/« 



where 



0J„ 



(36) 



This is nothing but the usual mass-radius relation for 
polytropes [43| extended to d dimensions [20], and writ- 
ten with the notations of biology. 

For n < ri3 there is one, and only one, steady state for 
each mass M and it is stable (global minimum of F[p] at 
fixed mass M). Its radius i?„ is determined by Eq. (|35|) . 
For 715 > n > 713 there is one, and only one, steady state 
for each mass M but it is unstable (saddle point of F[p] 
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at fixed mass). The index 77,3 is critical |20l |. For n = 713, 
steady state solutions exist for a unique value of the mass 



M c = S d 



"I «3/("3-l) 



A* 



(37) 



Their radius i?» is arbitrary and they are marginally sta- 
ble. This yields a family of density profiles of the form 



p(r)=p d ns (Zir/R*), 



(38) 

where po is related to i?* by Eq. (|3"T|). They all have 
the same mass (|3"T1) and it can be shown that their free 
energy is independent of i?» (see Appendix F of [20|). 
The invariant profile 0(f)" 3 is plotted in Fig. Q] In d = 3 
where 71,3 = 3, the critical mass is 



U, = 327TU.'., [ ^ 



3/2 



(39) 



It is found by solving numerically the Lane-Emden equa- 
tion that W3 = 2.01824... [l5(. Therefore, we obtain more 
quantitatively 



M, ■ : 202.8956... ( — 



.3/2 



(40) 



It can be convenient to introduce dimcnsionless variables 
or, equivalently, to take D = A = x = 1- I n that case, the 
only control parameter is the mass M. With this system 
of units, the critical mass in d = 3 is M c = 202.8956.... 




FIG. 1: Invariant density profile for the 3D generalized Keller- 
Segel model at the critical index n = 3 and critical mass M c . 

The critical mass (f3T)|) of bacterial populations is the 
counterpart of the Chandrasekhar limiting mass for white 
dwarf stars [llj. The analogy stems from the fact that, in 
the ultra-relativistic limit, white dwarf stars are equiva- 
lent to polytropes with the critical index n = 3 (in d = 3). 
Then, applying the theory of polytropes, Chandrasekhar 
[ll1 | obtains a unique value of the mass 



M, 



Chandra 



32tt 2 



1/2 



he 

c 1 77 



3/2 



1 



{pur 



(41) 



where h is the Planck constant, c is the velocity of light, 
G is the constant of gravity and H is the mass of the 
hydrogen atom (p is the molecular weight). In terms of 
the solar mass, it reads 



Mchar. 



= 0.196701. 



3/2 



1 



5.76M Q /7t 2 



{pHf 

In his more general treatment of par tially relativistic 
white dwarf stars, Chandrasekhar [44[ shows that the 
mass (|42[) represents a limit above which there is no equi- 
librium state [fpjsj. In that case, the system is expected to 
collapse leading ultimately to a neutron star or a black 
hole. 



V. THE TWO-DIMENSIONAL CASE 

In two dimensions (d = 2), the critical polytropic index 
713 — > +00. In that case, 7 = 1 and we recover the usual 
Keller-Segel model ©-(O corresponding to normal diffu- 
sion. The stationary solution is obtained by substituting 
the Boltzmann distribution (|10p in the Poisson equation 
([5|) . The resulting configurations have the same structure 
as isothermal stars in astrophysics. Their structure in d 
dimensions has been described in 0]. 

We consider spherically symmetric configurations. 
Defining 



Poe 



(43) 



where po is the central density, we find after simple alge- 
bra that tp is solution of the Emden equation 



1 d (td-iW 

e- 1 dt \* d$ 



(44) 



with ■0 = and ip' = at £ = 0. For d = 2, the density 
profile extends to infinity but the total mass is finite. The 



mass M = 2tt J + °° prdr is given by 



M = 



2ttD 



+ OO 



,-■<!> 



XA jo 

Using the Emden equation (|4"4l . we get 
2ttD 



M 



km £V/(0- 



(45) 



(46) 



In d = 2, the Emden function is known analytically 

1 



(1 + K 2 )' 

and we find that £-0' — > 4 for £ — > - 
unique value of the mass 



(47) 



-00. This leads to a 



8ttD 
ltA~' 



(48) 
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Therefore, unbounded steady states of the KS model in 
two dimensions exist for a unique value of the mass (|48[) 
and they are marginally stable ■ This yields a family 
of density profiles of the form 



p(r) = 



Pa 



(49) 



that are parametrized by the central density po- They 
all have the same mass (|48|l and it can be shown that 
their free energy is independent of po • The invariant 
profile e - ^ is plotted in Fig. [U 




FIG. 2: Invariant density profile for the two-dimensional 
Keller- Segel model at the critical mass M c . 



The critical mass of bacterial populations (|48|) is the 
counterpart of the critical mass or critical temperature 



M c .= 



4k R T 



knTc = 



GMr 



(50) 



Gm " 4 

of isothermal spheres in 2D gravity (see Ref. [l(| for a 
description of this analogy). The existence of a critical 
mass, or critical temperature, for systems described by 
the 2D Boltzmann-Poisson s yste m is known for a very 
long time in astrophysics 0,111, 51, S3, SI SI, HI, Eli, 
[52j and in the statistical mechanics of 2D point vortices 
[JIH, HHH^]. It has been rediscovered in the context of 
chemotaxis in [H, H, [H [H [H [H, [H, H| . Comparing 
Eq. (HI with Eq. ([33]) we find that for d = 2 and 
ii = 7T.3 — ► +oo, we have the limit 



lim n 3 w„ 3 



(51) 



This relation can also be obtained from Eq. (79) in [2Cl |. 
With this relation (f5Tj) . we can obtain the critical mass 
(|48f of the KS model in two dimensions as a particular 
case of the critical mass (|3"T1) of the GKS model when 
d — > 2. In the context of generalized thermodynam- 
ics, this corresponds to the limit 7 — ► 1 where we pass 
from the Tsallis (polytrope) distribution to the Boltz- 
mann (isothermal) distribution. In the present context, 
this corresponds to the passage from anomalous (7 > 1) 
to normal (7=1) diffusion. 



VI. CONCLUSION 

Our study tends to suggest that a notion of generalized 
thermodynamics can be useful in the context of chemo- 
taxis when the biological organisms experience anoma- 
lous diffusion. Therefore, the chemotactic aggregation 
of bacterial populations could be an important physical 
system where ideas of generalized thermodynamics ap- 
ply. The generalized Keller-Segel model (Ti"rJ)) - (rrT|) that we 
have introduced is related to the Tsallis form of entropic 
functional. This is a very rich model because it combines 
both elements of generalized thermodynamics and long- 
range interactions. Therefore, it enters in the general 
class of nonlinear mean field Fokker-Planck equations in- 
troduced by Chavanis [25| . In fact, we can develop an 
effective generalized formalism for an even larger class of 
entropic functionals than the Tsallis entropy. For exam- 
ple, we can consider a generalized stochastic process 



(52) 



where the diffusion coefficient is given by D(p) = 
xP{p)/p where P(p) is an almost arbitrary function 
which plays the role of an equation of state in the anal- 
ogy with barotropic stars in astrophysics [2{j. Then, we 
get the GKS model 



dp 
dt 



V • [ X (VP - pVc)] , 



Ac = -A/?. 



(53) 



(54) 



The equilibrium states correspond to the condition of 
hydrostatic equilibrium (TT4"|) . This system satisfies a for- 
malism of generalized thermodynamics [24| for a class of 
entropic functionals 



F = -- I pedv 



P 



(55) 



larger than the Tsallis entropic functional [59j. However, 
the Tsallis thermodynamics is convenient to describe de- 
viations from the Boltzmann thermodynamics in a simple 
manner. This leads to models that are still analytically 
tractable. For example, the GKS model (|16[) - (TTT1) can be 
studied in great detail [H, [H, [H[ . 

For the index n 3 , the GKS model (fTB|) - ([T7]) presents a 
critical dynamics involving a "universal" mass M c inde- 
pendent on the size of the system. In astrophysics, this 
result is well-known in d = 3 and is connected to the 
limiting mass of white dwarf stars discovered by Chan- 
drasekhar We can immediately transpose this re- 

sult to the biological context leading to the critical mass 
(1401) o f bacterial populations. This result was implicit 
in [20J and it has been developed in the present paper 
with emphasis. Thus, our paper reveals the analogy of 
the existence of a critical mass, found initially for white 
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dwarf stars, to exist also for bacterial populations driven 
by chemotaxis. Now, the question concerns the detailed 
description of the critical dynamics of the GKS model at 
n = 77,3 . This will be reported in a forthcoming paper [2^ | . 
From the present study, we anticipate that the critical dy- 
namics of the GKS model at 713 = 3 in d = 3 will be very 
similar to the critical dynamics of the ordinary KS model 
with 77,3 = +00 in d = 2 studied in [7J, |8|, |10( . This analogy 
will be confirmed in (29|. For M < M c , we find that the 
system evaporates (in an infinite domain) or tends to an 
equilibrium state (in a finite domain) corresponding to an 
incomplete polytrope confined by the box. For M > M c , 
we find that the system collapses. In a finite time, it 
forms a Dirac peak containing a mass M c surrounded by 
a halo that has a pseudo self-similar evolution. These 
results are similar to those found in d = 2. In conclusion, 



we can interpret the mass M c as a limiting mass above 
which the system undergoes chemotactic collapse. This 
strengthens the analogy with the Chandrasekhar mass of 
white dwarf stars. There is, however, a great conceptual 
difference between the two. Indeed, the Chandrasekhar 
mass is defined in terms of fundamental constants so it 
has a universal value. By contrast, the critical mass of 
bacterial populations depends on the parameters A, x 
and D that are not universal and that change from ex- 
periment to experiment. Furthermore, the index n = 3 
in d = 3 is special in astrophysics because it corresponds 
to the index of ultra-relativistic and completely degener- 
ate white dwarf stars ■ In biology, there is a priori no 
reason why the index n = n 3 should be selected in the 
dynamics of bacterial populations (except in d — 2 where 
it corresponds to a situation of ordinary diffusion). 
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[56] To be specific, we assume that the biological organisms 
experiencing chemotaxis are bacteria, but our results can 
be valid for other systems like amoebae, cells, social in- 
sects... The Keller-Segel model generically describes a 
wide variety of biological systems. 

[57] There are other methods to justify the GKS model (|f 6|l - 
<|17|) . It can be obtained from the master equation, assum- 
ing that the probabilities of transition explicitly depend 
on the occupation number (concentration) of the initial 
and arrival states [22, 40]. It can also be obtained in a 
strong friction limit of a generalized kinetic model taking 
into account inertial effects [241. l4lj. 

[58] The structure and the stability of (relativistic) white 



dwarf stars in d dimensions has been studied in 28]. It 
is found that the dimension d = 3 of our universe plays 
a very particular (marginal) role. 
[59] Note that the model (153ft - (1541 ) is still a particular case 
of the primitive Keller-Segel model CO)-©. However, the 
primitive Keller-Segel model does not satisfy a formalism 
of generalized thermodynamics. Apparently, the most 
general form of Keller-Segel model satisfying a formalism 
of effective generalized thermodynamics (effective tem- 
perature, Einstein's relation, //-theorem,...) is provided 
by Eqs. (134)-(135) of [p. 



